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ОПРЕДЕЛЕНИЕ ЭФФЕКТИВНЫХ СВОЙСТВ 
ПЬЕЗОКОМПОЗИТОВ НА ОСНОВЕ КОНЕЧНОЭЛЕМЕНТНОГО 
МОДЕЛИРОВАНИЯ В АСЕГАМ 


Разработана методика определения полного набора эффективных физических по- 
стоянных пьезокомпозитного материала нерегулярной структуры. Метод основан на 
решении набора динамических задач в конечноэлементном комплексе АСЕГАМ, для 
которого с этой целью разработан специальный модуль моделирования композици- 
онных материалов. Проведены расчеты по определению электрических и механиче- 
ских свойств пористой керамики. 

Ключевые слова: пьезокомпозиты, МКЭ, эффективные свойства, программные 
модули. 


Введение. Создание высокотехнологических пьезоэлектрических преоб- 
разователей энергии напрямую зависит от разработки новых пьезоактив- 
ных материалов. Один из путей получения этих материалов является раз- 
работка композиционных материалов с одной или несколькими активными 
фазами. Наиболее простой по составу среди пьезокомпозитов является по- 
ристая пьезокерамика с одной активной фазой. Различные способы изго- 
товления пористой пьезокерамики, которые позволяют варьировать харак- 
теристики материала, такие как объемное содержание, размер и распреде- 
ление пор, рассматриваются, например, в [1, 2]. Изучение электрических, 
акустических, упругих, пьезоэлектрических и диэлектрических свойств по- 
ристой пьезокерамики и выявление их зависимостей от микроструктуры 
пор способа изготовления и т.д. [3, 4] является актуальной задачей. Как 
показано в цикле работ Т.Г. Лупейко с соавторами [5-8], пористая пьезоке- 
рамика обладает рядом свойств, которые выгодно отличают ее от сплош- 
ной керамики и которые требуют своего объяснения. Причем эти свойства 
в значительной степени зависят от распределения фаз (содержание, 
открытость или замкнутость пор, размеры и форма пор и т.п.) 

Данная работа посвящена разработке методики расчета эффектив- 
ных свойств пористой керамики на основе решения ряда динамических за- 
дач для представительных объемов композитов и аналитических решений 
для тел с эффективными свойствами. При этом решение задач для пори- 
стой керамики проводится в специализированном конечно-элементном 
комплексе АСЕГАМ [9], в котором для этого разработаны модули моделиро- 
вания композитов нерегулярной структуры. 

Моделирование нерегулярных композитов в АСЕТАМ. На базе этого 
комплекса был создан модуль, позволяющий моделировать композитные 
материалы нерегулярной структуры, состоящие из нескольких фаз. Это мо- 
делирование происходит на уровне конечных элементов, которые наделя- 
ются свойствами определенной фазы, случайным образом (с равномерной 
функцией распределения) в соответствии с заданным объемным процент- 
ным содержанием фаз. Данный модуль встроен в решатель задач об уста- 
новившихся колебаниях составных упругих, электроупругих и акустических 
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тел, содержащих композиционные элементы. При этом предусмотрена воз- 
можность проводить множественные расчеты с определенным распределе- 
нием фаз. Для управления распределением фаз разработана специальная 
утилита в среде МсгооЁ Миа! Зи 2005 на языке С#, позволяющая ви- 
зуализировать это распределение. Пользователю предоставляется возмож- 
ность увеличивать/уменьшать масштаб изображения, а также включать/вы- 
ключать отображение номеров конечных элементов. Основной функцио- 
нальной особенностью данной утилиты является возможность изменения 
фазы конечного элемента при щелчке кнопкой мышки на его изображении. 
Также утилита позволяет отслеживать изменение объемного процентного 
содержания фаз. При моделировании пористой пьезокерамики предусмот- 
рена возможность определения свойства поры с помощью ослабления 
свойств основной активной фазы умножением на коэффициент понижения 
и превращением ее в пассивную фазу, а также возможна проверка связно- 
сти тела при рапределении пор. Первый вопрос, который возникает при 
моделировании пористого пьезокомпозита нерегулярной структуры, - это 
выбор достаточного коэффициента понижения и размера конечноэлемент- 
ной сетки, такого, при котором перераспределение с сохранением процент- 
ного содержания фаз не влияет на значения величин, найденных в ре- 
зультате расчета и используемых для определения эффективных свойств. 

Использование разработанного решателя с утилитой и выбор ука- 
занных выше параметров продемонстрируем на решении задачи о соб- 
ственных колебаниях в рамках линейной теории электроупругости [10] 
квадратной области (в системе координат Ох2 со стороной, равной / = 
0.009 м) в плоской постановке. Квадрат совершает колебания на первой 
толщинной моде в направлении вектора предварительной поляризации, 
параллельного оси О2. 

Краевые условия задачи имеют вид: 

при х=1//2 и=0,0,=0, Б,=0; 
при 2=+//2 0.=0,0,.=0,Ф=О, (1) 
где м,- компоненты вектора смещений; Ф - электрический потенциал; 
о’, - тензор напряжений; />, - компоненты вектора электрической 
индукции. 

В качестве материала сплошной керамики был выбран состав 
ЦТБС-3, свойства материала, моделирующего поры, задаются при помощи 
коэффициента понижения а, = а.10”, где < „ - упругие модули и плот- 
ность материала, задающего поры; @ - соответствующее механическое 
свойство материала сплошной керамики. 

Первая серия экспериментов была направлена на определение ко- 
эффициента понижения. На основании графика зависимости {[” относи- 


тельных резонансных частот первой толщинной моды ( = —^, У» 


о. 


частота первой толщинной моды пористого тела, Л. - частота первой 


толщинной моды сплошного тела) от степени коэффициента понижения И 
(рис.1) можно сделать заключение, что оптимальный коэффициент пони- 
жения составляет 107“. 


14 


Вестник ДГТУ, 2009. Т9. №1(40) 























038 083: 
095 0.821 
094 0.813 
Е 092 08- 
0.9 г ] 
088 В 
085 0.78] 
084 0.773 
0 05 1 15 2 25 3 35 4 45 5 3150 3200°3250°3300 3350 3400 3450 3500 
п м 


Рис.1. Зависимость относительной резонанс- Рис.2. Зависимость относительной резонанс- 
ной частоты от степени коэффициента по- ной частоты от количества узлов конечноэле- 
нижения ментной сетки 


Дальнейшие эксперименты были нацелены на определение наибо- 
лее оптимальной дискретизации тела. Исходя из графика зависимости от- 
носительных частот первой толщинной моды от числа узлов конечноэле- 
ментной сетки (рис.2) можно сделать вывод, что количество узлов сетки 
должно превышать 3300 узлов. Численные расчеты, проведенные при та- 
ком количестве узлов, показали, что перераспределение пор при сохране- 
нии их объемного содержания (8,66%) не влияет на первую толщинную 
моду 0,20 МГц (рис.З3, 4). 


















ГАД 
а 
Ч 
я 





СА 







ры 


ти 
Ч 






о 
ЧА 
(удУл 
А 
о 


БА 











№. 
о 
ы 
КУ 
Рух 
9 
АЯ 
р 
ты 
и 
и 


ты 
а 
м 
мы 
с 
ИЛ 
С 
к 
СУ 
к. 
точ 


- 

й 
ый 
ее 
р 


о 
ы 


ТАУ 
та 


| 
ДА 
р 
к 
че 


9 
ь 


ь. 
В 
Ри 





Г. 
ый 
7 
г 
й 
УЖ 


р 


АЛ 
ДУ. 
ь. 
и 
у, 
С 
УХ 
о 
© 
к] 
Е 
ы 
ы 
© 
ых 
ЕК 






Хх 
ы 
та 


> 
4 
та 
5 
А, 
С 
А. 
У 
РУ. 
У 
й 


- 


а 74 ро 
АА 
о СА 
а, 
Я 


Г 
ы 
У 
м 
т 
КУ». 
ь 















та ть а ть 
Е 






"ЗКУАУА. А. АА 
ИУ ЗЕ 








Рис.3. Первоначальное распределение и перераспределение пор в 8,66%-ном композите 
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Рис.4. Первая толщинная мода 
Определение эффективных свойств с помощью динамических за- 
дач. Разработанный решатель применен в задаче определения эффектив- 
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ных упругих и пьезоэлектрических свойств композитной пьезокерамики не- 
регулярной структуры. Одним из способов описания композитных пьезома- 
териалов является замена их эквивалентными однородными составами с 
некоторыми эффективными характеристиками. Задаче определения полно- 
го набора материальных констант пьезокерамики уделяется большое вни- 
мание в литературе [11, 12]. Несмотря на ряд разработанных методик 
[13, 14], определение этих свойств и выявление их зависимостей от рабо- 
чих параметров (частоты колебаний, температуры и т.д.) является актуаль- 
ной задачей. 

Далее рассматривается способ определения эффективных упругих 
и пьезоэлектрических свойств композитной пьезокерамики нерегулярной 
структуры, основанный на динамической эквивалентности ее однородной 
пьезокерамике. Для этого рассматриваются собственные колебания некото- 
рых представительных объемов композита и находятся их резонансные и 
антирезонансные частоты. Их вычисления проводятся в конечноэлемент- 
ном комплексе АСЕГАМ. Эти частоты для ряда одномерных и двумерных 
движений, описывающихся аналитически, являются исходной информаци- 
ей для определения эффективных констант. При этом предполагается, что 
ряд характеристик, таких как диэлектрические проницаемости и плотности, 
могут быть измерены по известным экспериментальным методикам для од- 
нородных материалов [11] (в численном эксперименте они вычислялись по 
соотношениям теории смесей). В качестве таких мод колебаний рассматри- 
ваются: 

А) одномерные поршневые движения, параллельные вектору пред- 
варительной поляризации (толщинные моды); 

Б) колебания перпендикулярно этому вектору протяженных тел со 
свободными границами и поршневые движения в этом направлении (про- 
дольные моды); 

В) сдвиговые колебания в направлении вектора поляризации. 

В рамках линейной теории электроупругости [10] рассматриваются 
собственные колебания прямоугольных тел в плоской постановке. В каче- 
стве материала сплошной керамики, как и ранее, был выбран ЦТБС-3, име- 
ющий следующие механические свойства: 








15.1 7.9 80 0 0 0 
15.1 8.0 ооо 
С- 13.6 0 оо 
29 0 0 
29 0 
симм 3.6 
0 0 0 0 154 0 
Е= 0 0 0 154 0 0; 
- 7.9 - 7.9 17 0 оо 
14288.5 0 0 
ё = 0 14288.5 о, 
0 0 11328 
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где С. 10° Н/м? - матрица упругих модулей; Н - матрица пьезомодулей; 


с . 109 12 Ф/м - матрица диэлектрических проницаемостей. 
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Решение задачи с граничными условиями (1) ищем в виде: 
и, (х,2)= 0, и, (х,2)= И.(2), ф(х,=) = Ф(2). (2) 
Уравнения, определяющие частоты резонанса и антирезонанса име- 

ют вид: 
1.2 
[5 771 °.1 — ь о, ия а (3) 
2 
где. т” = ай К, = би р= 13,5, @,, @ 
ё ззСзз р 





- собственные частоты 





а 


резонанса и антирезонанса в толщинных колебаниях, =,, ©, с, - 
эффективные компоненты тензоров диэлектрической проницаемости, 
пьезоконстант, упругих постоянных. 
Краевые условия задачи Б, имеют вид: 
при х= +//2 0, =0,06,=0, р =0; 
при ==+й/2 0.=0,0.=0,Фф=0. (4) 
В предположении, что />> Я иод..(х,2) = О решение задачи 
ищется в виде: 
и, (х,2)= И, (%); и, (х,2) = И. (2); Фе, 2) = Ф(=). (5) 
Уравнения, определяющие частоты резонанса и антирезонанса, 








имеют вид: 
К 
9, = (2п- 1 РЕ 2; 
1 1 1 1 (1- Р)а- =^)с 
2 в, —1- 8° |= © 1 - 2? 
+ м2 р о = в. (5) 
2 
ДЕ 2? = С ‚е= е с с= —. 
С Сзз а. би 


Краевые условия задачи Б2 имеют вид: 
при х= +//2 0, =0,0,=0, Б=0; 


при 2=+й/2 и=0,0,=0,ф=0. (7) 
Решение задачи ищется в виде: 
и, (%,2) = И, (%); и, (х,2)=0; Ф(х,Р)=0. (8) 
Уравнения, определяющие частоты резонанса, имеют вид 
К 
@, = (2п- Л г (9) 
Краевые условия задачи В имеют вид: 
при х=1//2 0, =0,0,=0,ф=О0; 
при 2=+й/2 06.=0,06,=0, О, =0. (10) 


Решение задачи ищется в виде: 
и, (х,2)= 0; и, (х,2) = Ц .(х); Ф(х,2) = Ф(Х). 
(11) 
Уравнения, определяющие частоты резонанса и антирезонанса име- 
ют вид [15]: 
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2-0. 7=е, (12) 


В № | 
о, 51 —®, аЕ Эа: 
2 Ара А в 
Таким образом, система уравнений (3), (6),(9), (12) на каждой гармо- 
нике представляет собой систему семи уравнений относительно семи неиз- 


вестных: С\1› Сз› Сзз› С55› е1› е.з› е,5 - эффективных констант плос- 


кой задачи теории электроупругости для материалов класса бтит. При 
этом предполагается, что диэлектрические проницаемости измерены в экс- 
перименте. Для их численно аналитического нахождения могут быть 
рассмотрены аналогичные задачи для осесимметричных представительных 
объемов. 
В качестве результатов численного эксперимента приведем данные 
о нахождении модуля упругости с.- и пьезоконстанты е.. на основе 
уравнений (3) и результатов расчета частот резонанса и антирезонанса 
представительного объема 10,0% и 23,0% пористой керамики 
нерегулярной структуры в комплексе АСЕТАМ (таблица). 
Результаты численного эксперимента 





где 7 = 


Константы Эффективные константы пористой керамики 
пьезокерамики ЦТБС- Объемное содержание Объемное содержание 


. пор-10% _Х.. =202.2 пор-23.0% _Х’. =165.5 


кГц, _Д., =220.4 кГц кГц, Д., =184.8 кГц 


е.. 17, 7Кл/ м” 13,717Кл/ м” 11,177 Кл / м” 
с 13,6х10° Н/м? 8,131х10° Н/м” | 4,713х 10° Н/м? 





Заключение. Разработанная методика позволяет определять эффектив- 
ные свойства пористой пьезокерамики для известного состава пьезокера- 
мики и процентного содержания пор на основе проведения численных рас- 
четов в конечноэлементных пакетах, в которых возможно моделирование 
композитов нерегулярной структуры. Проведение натурных экспериментов 
для макрообъемов пьезокомпозитов по известным методикам [11] в тех 
случаях, когда это возможно осуществить, позволит уточнить данные чис- 
ленного эксперимента. 

Работа выполнена при частичной финансовой поддержке РФФИ 
(коды проектов 07-08-13589, 07-08-12193). 
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